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Abstract 

This  paper  presents  comparison  between  full  order  and 
minimum  order  observer  controller  design  using 
MATLAB  environment  As  a matter  of  fact,  it  has  shown 
how  to  implement  a dc  motor  state  space  model  and  its 
observer.  Different  observer  gain(s)  are  determined  by 
selecting  different  Eigen  values  for  the  observers.  Both 
designed  controllers  are  tested  using  root  locus  method 
by  varying  different  observer  poles.  Having  full 
understanding  of  the  observer  controller 
implementation(s),  students  and  engineers  will  fill 
confident  to  use  this  observers  and  observer  based 
controllers  in  numerous  engineering  and  scientific 
applications. 
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1.  Introduction 

In  this  paper  we  shall  consider  a problem  of  designing  a 
regulator  systems  by  using  the  the  pole-placement-with  - 
observer  approach.  We  shall  use  the  following  design 
procedure. 

1 . Derive  a state  space  model  of  a plant.  In  this  paper  a DC 
motor  is  considered  as  a plant  and  armature  current  and 
speed  of  the  dc  motor  are  taken  as  state  variables. 

2.  Chose  the  desired  closed-loop  poles  for  pole  placement 
and  chose  the  desired  observer  poles. 

3.  Determine  the  state  feedback  gain  matrix  K and  the 
observer  gain  matrix  Ke. 

4.  Using  gain  matrix  K and  Ke  obtained  in  step  3,  derive 
the  transfer  function  of  the  observer  controller.  Root  Locus 
methd  is  used  to  analyse  the  stability  of  the  controller 
based  on  location  of  the  observer  poles. 

Observers  are  used  these  days  not  only  for  the  purpose  of 
feedback  control,  but  also  in  their  own  right  to  observe 
state  variables  of  a dynamic  system,  which  can  be  an 
experiment  in  progress  whose  state  has  to  be  monitored  at 
all  times.  However,  in  the  most  recent  editions  of  several 
standard  undergraduate  control  system  textbooks  we  can 


find  the  coverage  of  full-order  and  even  reduced-order 
observers  [2] -[4]. 


2.  DC  Motor  Modeling  Using  State  Space 
Analysis 


The  different  equations  related  to  DC  motor  are  given 
below  [11], 
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Where  ec(t)  = armature  voltage,  em  (t)  = back  emf,  L(t)  = 
armature  current,  Tit)  = developed  torque,  9(f)  = motor 

■iS'Ti 

shaft  angle,  . = co(t)  = shaft  speed,  J = moment  of 

inertia  of  the  rotor,  B = viscous  frictional  constant,  Lm  = 
inductance  of  armature  windings,  R m = armature  winding 
resistance,  A^r=  motor  torque  constant,  Km  = motor 
constant. 


Here  the  motor  speed  co(t)  is  controlled  by  varying  the 
armature  voltage  ^(t).  Hence  ec(t)  is  the  input  variable 
and  co(t)  is  the  output  variable. 

i s ' n 

We  chose  as  the  state  variables  x^f)  = oo(tj  = and 


x2  (tj  = !'<z  (0  (5) 

The  state  equations  will  now  be  derived  by  using  above 
equations. 
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Hence  state  model  of  dc  motor  is  derived  from  equations 
(6),  (7)  and  (8)  as  follows 
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3.  DC  Motor  State  Model  Using  Motor 
Parameters 


Let,  the  motor  parameters  (coefficient  of  differential 
equations)  are  assigned  to  be  = 0.5  H,  Kt  = 0.01  N- 
m/A,  Km=  0.01  V-sec/rad,  J = 0.01  kg-m2,  B = 0.1  N-m- 
sec/rad,  Rm=  1 Q. 

Thus  the  state  model  of  dc  motor  is  derived  using  motor 
parameters  and  equation  (9)  and  (10)  as  follows: 


:;ri 


(12) 

4.  Full  order  observer  design 

The  theory  of  observers  originated  in  the  work  of 
Luenberger  in  the  middle  of 1960s,  [7] -[9]  .According  to 
Luenberger,  any  system  driven  by  the  output  of  the  given 
system  can  serve  as  an  observer  for  that  system.  Consider 
a linear  dynamic  system  with  unknown  initial  value  of  x 

x = Ax  -f  Bu  , y = Cx  (13) 

We  define  the  mathematical  model  of  the  observer  to  be 
£ = (A-KeC)X+  Bu+  Key  (14) 

Where  X is  the  estimated  state  and  C X is  the  estimated 
output.  The  nxl  matrix  Ke  is  called  state  observer  gain 
matrix.  Hence  the  observer  error  equation  is  defined  by 

x - x = (A  - KeC)  (X  -x) 

e = (A-KeC)  e (15) 


Where  (X  -x  ) = e = observer  estimation  error  vector.  Thus 
the  dynamic  behavior  of  the  error  vector  depends  upon  the 
Eigen  values  of  A - KeC  . 

If  the  observer  gain  Ke  is  chosen  such  that  the  feedback 
matrix  A-KeC  is  asymptotically  stable  (has  all  eigen 
values  with  negative  real  parts)  then  the  estimation  error  e 
will  decay  to  zero  for  any  initial  condition  of  e.  This 
stabilization  requirement  can  be  achieved  if  pair  (A,C)  is 
observable. 


The  system  (14)  under  the  perfect  state  feedback  control, 
that  is 
u = —Kx 

has  the  closed  loop  form 

i = (A  — BK)x  (16) 

So  the  eigen  values  of  the  matrix  ,4  - B K are  the  closed 
loop  eigen  values  under  perfect  state  feedback. 

For  state  feedback  control  r , It  is  known  that 
u = -Kx 

With  this  control,  the  state  equation  (13)  becomes 
k = Ax-  BKx 

= U - + 

= Gi  — BK}x  +BKe  (17) 

Where  X -x  = & = observer  estimation  error 
Combining  equations  (17)  and  (15),  we  obtain 


x]_\A  - BK  BK  1 m 

b\l  0 iJ-ALcJl&l 


(IB) 


The  above  equation  describes  the  dynamic  of  the  observed 
state  feedback  control. 


Ackerman’s  formula  is  used  to  write  the  program  in 
MATLAB  command  window.  The  programs  are  given 
below  [12]. 

» %State  feedback  matrix  K design  using  pole- 
placement  technique 

A = [-10  l;-0.02  -2]; 

B = [0;  2];C  = [10];D  = [0]; 

J = [-2+j*l  -2-j * 1 ] ; 

K = acker  (A,B,J) 

K = 


32.4900  -4.0000 


» % State  observer  gain  matrix  Ke  design 


» A = [-10  l;-0.02  -2]; 

B = [0;2] ; C = [1  0];  D = [0]; 
L=  [-9-10]; 
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Ke  = acker  (A',  C\  L)' 

Ke  = 

7.0000 

55.9800 

5.  Full  Order  Observer  Based  Controller 
Design  for  DC  Motor 

The  transfer  function  of  the  observer  controller  is  given  as 

$t-7iir-'*s,-A-K'c+ BK>~,K' 

(19) 

The  above  transfer  function  can  be  represented  by 
following  block  diagram  model. 


System 


Figl:  Block  diagram  of  system  with  observed  state 
feedback 

The  following  MATLAB  program  is  used  to  produce  the 
transfer  function  of  the  full  order  observer  controller  [12]. 

» A=[-10  l;-0.02  -2]; 

B=[0;2];C=[1  0];D=[0]; 

K=[32.5  -4]; 

Ke=[7;56]; 

AA=A-Ke*C-B  *K; 

BB=Ke; 

CC=K; 

DD=0; 

[num,  den]=ss2tf  (AA,BB,CC,DD) 


0 3.5000  35.5600 

den  = 

1.0000  11.0000  19.0200 

Hence  the  transfer  function  of  the  full  order  observer 
controller  is 

3.5s  -f  35.56 

s2 -h  11s  -h  19.02  (20) 

6.  Full  Order  Observer  controller  analysis 
based  on  Root  Locus 


Root  Locus 


Fig2:  The  root  locus  plot  of  the  full  order  observer 
controller  with  observer  poles  at  s =-9  and  s =-10. 


Root  Locus 


Fig3:  The  root  locus  plot  of  the  full  order  observer 
controller  with  observer  poles  at  s =-14  and  s =-15 


num  = 
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It  is  seen  from  the  above  two  figs  that  If  we  place  the 
observer  poles  far  to  the  left  of  the  jw  axis,  the  observer 
controller  become  unstable. 

7.  Minimum  Order  Observer  Design 

Consider  equation  (14)  where  the  state  vector  X can  be 
partitioned  into  two  parts  x^  (a  scalar)  and  xb  (a  vector). 
Here  the  state  variable  xG  is  equal  to  the  output  y and  thus 
can  be  directly  measured  and  x^  is  the  unmeasurable 
portion  of  the  state  vector.  Then  partitioned  state  and 
output  equations  become 


y=[i  010  (22) 


From  equation  (15)  the  equation  for  the  measured  portion 
and  unmeasured  portion  of  the  state  becomes. 


+ 3,U 

(23) 

xs  + Aibxb  +3ju 

(24) 

Here  equation  (18)  and  (17)  are  known  as  ‘state  equation’ 
and  ‘output  equation’  for  the  minimum  order  observer. 

The  equation  for  the  full  order  observer  is 

x = (A-KeC)X+  Bu+  Key  (25) 

Then  making  the  substitution  of  table  1 into  last  equation 
we  obtain 

*b  = ( Abb  ~ KA<zb)Xb  + ^£,flX£  + B& u + ks  Aeb  ** 

(26) 


Full  order  state 

Minimum  order  state 

observer 

observer 

X 

A 

A 

^ bb 

Bu 

u 

y 

xa  - 4::x:  - B.-u 

c 

kf  (nX  1 matrix) 

k€  [(n— 1)  Xl  matrix] 

Table4:  List  of  necessary  substitutions  for  writing  the 
observer  equation  for  the  minimum  order  state  observer 


~ — (Afeb  - “ Ab  ) (27) 

Define  e = x-L1  - x% 

Then  equation  (21)  becomes 

* = (Abb  - Mab  ) e (28) 


This  is  the  error  equation  for  the  minimum  order  observer. 
The  dynamic  behavior  of  error  vector  depends  upon  the 
Eigen  values  of  Abb  - . 

Now  the  characteristic  equation  for  the  minimum  order 
observer  is  obtained  the  equation  as  follows. 

1 5/  - Abb  -b  k6  4-bl  = (s  -mL)  (s -mz)...  (s  -mR_ i)  = 
0 

Where  ,m3 , mn_  L are  desired  Eigen 

values  for  the  minimum  order  observer.  Suppose  the 
desired  location  of  the  Eigen  value  for  the  minimum 
observer  is  at  s = -9. 

The  minimum  observer  gain  matrix  is  designed  by  the 
following  MATLAB  program. 

Ackerman’s  formula  is  used  to  write  the  program  in 
command  window.  The  program  is  given  below  [13]. 

» % State  observer  gain  matrix  Ke  design. 

Aab  = [1];  Abb  = [-2]; 

LL  = [-9]; 

» Ke  = acker  (Abb',  Aab',  LL)' 

Ke  = 


7 


8.  Minimum  Order  Observer  based 
Controller  Design  for  DC  Motor 


The  transfer  function  of  minimum  order  observer 
controller  is  given  as 


L-i  ;-r-  _ nurr. 


t (5/  - +5] 


(29) 


The  above  transfer  function  can  be  represented  by 
following  block  diagram  model. 


By  subtracting  equation  (26)  from  equation  (24)  we  obtain 
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Transformation 


Fig5:  System  with  observed  state  feedback,  where  the 
observer  is  minimum-order  observer 


Hence  the  transfer  function  of  the  minimum  order  observer 
4.5  (£+2.8) 

controller  is  — : (30) 

(5+1)  v ' 

9.  Minimum  Order  Observer  controller 
analysis  based  on  Root  Locus 


Root  Locus 


Fig6:  The  root  locus  plot  of  the  minimum  order 
observer  controller  with  observer  pole  at  s =-9 


Where  A = A - F kb 

S=B-F(ka+kbks) 

C = -kb 

S = - (ka  +kbka) 

F =B-b  -k6B.- 

The  following  MATLAB  program  is  used  to  produce  the 
transfer  function  of  the  minimum  observer  controller  [13] 
A=[-10  1;-0.02  -2]; 

B = [0;2];C  = [10];D=  [0]; 

Aab  = [1];  Abb=  [-2]; 

Ke  = 7;  K=[32.5  -4];  Kb=-4; 

Ka=32.5; 

Aaa  = -10;  Aba  = -0.02; 

Ba  = 0;  Bb=2; 

Fhat  = Bb-Ke*Ba; 

Ahat  = Abb-Ke*Aab; 

Bhat  = Ahat*Ke+Aba-Ke*Aaa; 

Atilde  = Ahat-Fhat*Kb; 

Btilde  = Bhat-Fhat*(Ka+Kb*Ke); 

Ctilde=  -Kb; 

Dtilde  = - (Ka+Kb*Ke); 

[num,  den]  = ss2tf  (Atilde,  Btilde, -Ctilde, -Dtilde) 
num  = 

4.5000  12.5800 


den  = 

1 1 


Root  Locus 
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Fig7:  The  root  locus  plot  of  the  minimum  order 
observer  controller  with  observer  pole  at  s =-14 

It  is  seen  from  the  above  two  figs  that  If  we  place  the 
observer  poles  far  to  the  left  of  the  jw  axis,  the  observer 
controller  become  unstable 

10.  Conclusion 

This  paper  has  shown  in  detail  how  to  implement  full- 
order  and  minimum-order  observer  controllers  based  on 
root  locus  in  MATLAB  environment  and  presented 
corresponding  fundamental  derivation  and  results  with  the 
help  of  dc  motor  state  model  and  the  conclusion  is  if  the 
observer  controller  becomes  unstable,  move  the  observer 
poles  to  the  right  in  the  left-half  s plane  until  the  observer 
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controller  becomes  stable,  Also  the  desired  closed-loop 
pole  location  may  need  to  be  modified. 
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